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Abstract. In this article we rederive the Boost operators in Coulomb-Gauge Yang-Mills theory employing 
the path-integral formalism and write down the complete operators for QCD. We immediately apply them 
to note that what are usually called the pion square, quartic... charge radii, defined from derivatives of 
the pion form factor at zero squared momentum transfer, are completely blurred out by relativistic and 
interaction corrections, so that it is not clear at all how to interpret these quantities in terms of the pion 
charge distribution. The form factor therefore measures matrix elements of powers of the QCD boost and 
Moeller operators, weighted by the charge density in the target's rest frame. In addition we remark that the 
decomposition of the 77' wavefunction in quarkonium, gluonium, ... components attempted by the KLOE 
collaboration combining data from 4> radiative decays, requires corrections due to the velocity of the final 
state meson recoiling against a photon. This will be especially important if such decompositions are to be 
attempted with data from J/t/j decays. 



PACS. ll.30.Cp Lorentz and Poincare invariance 
13.40.Gp Electromagnetic form factors 



1 Introduction 

1.1 The pion form factor 

Form factors traditionally encode the structure of a com- 
posite target as accessible from elastic reactions. For a 
scalar target, we can denote the elastic form factor simply 



as F(Q2), where 



p — p'. It is conventionally 



normalized as F(0) = Q, the particle's charge. 

For a non-relativistic particle with unit charge (Q = 1, 
~ 0), an expansion around zero momentum transfer 
allows for a physical interpretation of the form factor in 
terms of its rest frame charge density p(r), given by 
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F(Q^) = 1 - ^(r').Q' + 1:t(OpQ' + O^. (1) 
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Here we used the notation of Ref. [T]. In hadron physics 
one often quotes also the curvature of the form factor [5] 
via 



F(Q^) = l--(r- 



(2) 



Comparing with Eq. ([T|) we see that "^^y = 5t(i'^)- 

Beyond the charge normalization / d^rp(r) = (l)p = 1, 
the derivative at the origin provides, in non-relativistic 
quantum mechanics, the target charge radius 

-6 ^ . (3) 



dQ2 



3=0 



13.20.Gd Decays of J/ip, and other quarkonia 



It is known from nuclear physics, but sometimes ignored 
in the hadron physics community, that the non-relativistic 
interpretation should be modified as effects of boosting 
the pion wave function begin to appear. Since the pion is 
such a light target, 1/m^ corrections are even larger than 
the actual (r^) measured experimentally or computed on 
a lattice. 

The situation is even worse for the curvature, that has 
been the object of our recent focuse. Using both chiral per- 
turbation theory and dispersion relations, we have found a 
reliable value for the curvature of 4.0 ± 0.5 GeV~^. Many 
other authors quoted in [3] have also dedicated time to ex- 
tracting this curvature. It behooves us to examine whether 
this number -and the more common (r^), have any inter- 
pretation in terms of the target wavefunctions. We show, 
indeed, that this is far from trivial. In Fig. [1] we replot the 
form factor data from the NA7 experiment at CERN [i], 
omitting the few points at higher momentum transfer (as 
they have large error bars that might still be underesti- 
mated). One can then obtain a direct fit that is sensitive to 
the "square radius" and curvature. From the parabolic fit, 
one can derive F'(0) = -2.017±0.047 GeV-^, or in units 
of square radius, (r^) = 0.48 ± 0.01 fm^. The curvature is 
less well determined, ~ 8 ± 1 GeV~^. Adding a small 
systematic error to the data the uncertainty in the curva- 
ture could be as high as 2 GeV~^, still somewhat higher 
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Fig. 1. Linear and parabolic fit to the elastic pion form fac- 
tor at small momentum transfer. The linear fit has x^/dof ~ 
77/32. The parabolic fit has xV^of ~ 22/32. 

than the result obtained from the Omnes representation 
quoted above. 

Several works in atomic and nuclear physics have pointed 
out the necessity of 1/m corrections when talking about 
the charge radius of a composite target as extracted from 
electron scattering or from the isotope shift [5] , for exam- 
ple the so called Darwin- Foldy terms jB] . In general much 
of the literature features idle discussion about where the 
various effects are to be adscribcd, to the radius, or to its 
extraction from electron scattering, or to final state rescat- 
tering, etc. In the end, a lot of the discussion is linguistic. 
In nuclear physics one is not quite yet in possession of a 



full Hamiltonian for the fcw-nucleon problem, so the cor- 
rections are worked out on a semiempirical basis [7] and 
in atomic physics the corrections due to both relativity 
and interactions are small, and the first order ones can be 
accessed by means of simple quantum-mechanical reason- 
ing, although some works exist providing first-principles 
computations [8]. 

In hadron physics however, the situation is radically 
different. Corrections are large, since energies involved 
match particle masses (the most salient example being the 
pion) and interactions strong. In addition one is in posses- 
sion of the Hamiltonian of the fundamental theory, QCD, 
and in this article we seek the first brick in the direction 
of gaining understanding of form factors in terms of the 
canonical quantization of the theory in Coulomb gauge. 

Much is known already in light-front pj and point- 
form quantization pilllllj . with recent insight from AdS- 
QCD [H]. However a good discussion of boost operators 
of QCD in a phenomenological context is nowhere to be 
found in the literature for canonical. Coulomb gauge quan- 
tization. In related quark models there has been indeed 
interesting work within the Bakamjian- Thomas construc- 
tion |13j. Here we attempt to continue stimulating the 
discussion and bring it a step closer to QCD. 

1.2 0-radiative decays and Fock space expansion 

The KLOE collaboration [Tl] has proposed to employ their 
excellent data set of radiative decays to pseudoscalar 
mesons, ^ — >■ 777, 777', to measure the gluonium con- 
tent of the T]' as well as the pseudoscalar mixing angle. 
The theoretical assumptions behind the analysis are sim- 
ple considerations about the flavor structure of the cou- 
plings of gluonium, quarkonium, etc. 

While there is no objection to the model analysis of the 
pseudoscalar mixing angle [12], we believe that a point 
has been missed by the community in what regards the 
gluonium content of the rj' . This is the fact that the Fock- 
space expansion for either of the 77 or rj' 

\v) =ci|qq)+c2|00)+c3|qqg)-f ... (4) 

is dependent on the reference frame. When one makes an 
assumption about the gluonium (gg and larger number of 
gluons) content of a meson, one is presumably referring to 
the rest-frame of the meson, since the concept of quantum 
mechanical state requires a quantization surface that is 
presumably taken to be t = in the rest frame. But, 
although the rj' meson is produced non-relativistically in 
(j> radiative decays, the analysis hinges also on decays to 
the T] meson, and this is produced with a velocity v ~ 0.55. 

As we will show, even if there was no gluonium at all 
in either of the two mesons, one should expect a gluonium 
content of the 77 meson boosted to such reference frame 
of order sin(0p)?;^/2 ~ 0.1, given that the pseudoscalar 
mixing angle seems to be close to 38 degrees. This is the 
well-known effect that the boost operators involve inter- 
actions and change the particle content. 

This is not necessarily a fatal flaw of the analysis, but 
a call of attention that, while effective hadron Lagrangians 
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are explicitly Lorentz-invariant, such Fock-space decom- 
positions are tied to a reference frame where canonical 
QCD quantization is carried out, and this needs to be 
specified and consistently handled. 

The necessity of specifying a reference frame will be- 
come more accute as the natural next step in the analy- 
sis [l6lfT7] will be to employ the data basis for radiative 
J /ip decays accumulated by BES and others, and given 
the mass of 3097 MeV of the cliarmonium ground state, 
the velocity of the produced mesons will now be decidedly 
relativistic. For J/-0 —J- 77?, = 0.94, and for J /ip — 777', 
Vrj' = 0.83, that are indeed significant. 

The rest of this paper is organized as follows. In sec- 
tion [2] we motivate the necessity of constructing the Boost 
operators with a careful setup of the form factor for a 
scalar target. Then in subsection 12.11 we particularize to 
the case of the pion, the lightest hadron where relativis- 
tic corrections are most significant. Subsection 12. 21 is then 
dedicated to dispelling some misconception about the charge 
radius interpreted in the Breit frame. 

Although pieces thereof can be found in the literature, 
it is convenient to write down the boost operators of QCD 
in Coulomb gauge for future reference, which we do in sec- 
tion [3] There we also present a relatively easy derivation 
based on modern functional methods in subsections l3.1l3."^ 
and 13.31 

Finally, section|4]presents the application of these boost 
operators to 777 (or similar decays) where it has 
been until now ignored, and section [5] presents our final 
remarks. 



2 Rigorous interpretation of the form factor 
in canonically quantized field theory 

The form factor can be expressed, in non-relativistic nor- 
malization for the charge density, as 



(P'lf|p>in = 



(P + P')^ 



F(q2). 



(5) 



The energy-momentum p^ = (m7r,0) corresponds to a 
pion at rest and time t = —00 before the momentum 
transfer by the virtual photon, q^, boosts it to a frame 
with p'. 

The charge density is taken at t = 0, that we employ 
as quantization surface to define the equal-time commuta- 
tion and anticommutation relations for field operators. To 
undertake any computation one needs to propagate the 
initial state |p)in to t = 0. The propagators that accom- 
plish this operation are usually called Moellcr operators in 
the interaction picture, = U(t = —00, 0). Likewise the 
propagation of the final state to the quantization surface 
is carried out by f2\_ = U'f(t = 0, 00). The product of both 
operators would reconstruct the S matrix, !7]_i7_ = S, but 
here an insertion of the current operator occurs at time 
t = 0, see Eq. ^ below. 

The final-state pion is boosted, and to be able to use 
its wavefunction in the rest frame (where the charge den- 
sity is defined) one needs to employ the boost K operator 



from QCD 

i„(p'l =i„ (Ple"*'*-^ (6) 
that we construct in Sec. [3] below. In terms of these oper- 
ators the form factor can be expressed via 

(P + P')'' 



i„(p|e-'K-«4.7^(0)f2_|p)in = 



-F(q^) (7) 



(the MccUer operators introduce interactions both sides of 
the current insertion). An alternative form of i7^j^(0)51_ 
is for example T (^j^(0)c*/'i*''~i) . 

The interpretation of the form factor in non-relativistic 
quantum mechanics (where it comes about in the Born 
approximation by expanding the Fourier transform of the 
potential causing the scattering) 



FNonRci(|qn = Q-^|qP(r2) 



1 



q|'(r') 



(8) 



3!'"' ^ ' 5!' 
in terms of charge radii, is only a limiting case that can 
be recovered from Eq. ([7|) and we will do it shortly. 

Let us now obtain an equivalent q^/mj expansion of 
Eq. ([7|) in powers of the momentum transfer to match 
Eq. dH]). The rapidity parameter in Eq. © depends on the 
transferred momentum through 



V = tanhC =^ C = ^ log 
From the above equation we can approximate 



C ~ w H h . . . 0{v^) with 



3 q^ 
4n4 



(9) 



(10) 



Expanding also E' in terms of q, and making q oc 
(hz = kz/|kz|) infinitesimal so that one can truncate the 
Taylor expansion of the boost 

exp (^iK^C^^ 



4! 



we obtain 



F(q2 



(11) 




(12) 



4j"f7_|p) 



(terms odd in are absent for a scalar target because of 
z — >■ — z refiection symmetry). We then have 

F(q2) ^ F(0) + q2F'(0) -I- ^qV(O) + . . . (13) 

where the derivatives arc now 

F(0) = (p|r!^jOr!_|p), 

F'(0) = 2SjF(0) + ^(p|K2f];jOi7_|p) (14) 



4 



Gomez Rocha, Llanes-Estrada, Schiitte and Vill 



2.1 No square radius interpretation for the pion 

To reproduce Eq. ([5]) , one first invokes the impulse approx- 
imation, by which the quarks from the j'' current density 
directly enter the state's wavefunction, in effect neglecting 
the Moeller operators (rescattering) . This is known from 
nuclear physics to be a poor approximation, and many 
works in hadron physics have also lifted it [T8| since it 
fails, for example, to reproduce simple instances of Vector 
Meson Dominance. 

In a second step, one takes the non-relativistic limit 
by sending the mass denominators m^ — >■ oo, and substi- 
tutes the Lorentz boost by its Galilean equivalent. Since 
a Galilean change of reference frame is generated by 



U[v] = cxp (— ini^f • v) 



(15) 



the relevant boost operator in quantum mechanics be- 
comes K — >■ Kcaiiican = —iTiTri- Upgrading it to non- 
relativistic field theory for quarks, and summing over spin, 
color and flavor, the boost operator is then 



K 



Galilean 



E 



(16) 



(In terms of quadrispinors including antiparticles, there 
would be an additional Dirac beta matrix (— imq. )/39/i9k' 
in the operator). 

In this case, one obtains for the form factor a non- 
relativistic expression in terms of the wavefunction of the 
target. If the target |H) is taken as a two-body scalar with 
non-relativistic wavefunction =^(r), and two particles of 
mass mi, m2 in the center of mass frame, with relative 
position r, and opposite charge (as corresponds for exam- 
ple to the non-relativistic Hydrogen atom or the neutral 
ttq in the quark model) 



; J dW|^(r)|= 



mi - m2 
mi + m2 



(17) 



However, returning to Eq. ([T^ for the case of the 
pion, l/2m^ ~ 26 GeV~^ is not negligible in any sense 
against — F'(0) ~ 2 GcV~^ taken from the experimental 
data quoted above. There need to be large cancellations 
among the two terms in the second of Eq. (fH)) . Its sec- 
ond non-constant term, that in the non-relativistic limit 
would correspond to the square radius, is therefore of or- 
der 20 GeV~^. One needs to conclude that, in the case 
of the pion, there is no connection whatsoever between 
the derivative of the spacelike form factor and the square 
radius of the rest frame charge distribution. 



2.2 Breit frame 

In the study of the proton form factor, the Breit frame 
became very popular |191l20j . One reason is that, for a spin 
1/2 target, a separation of the charge and magnetization 
densities is possible by employing the Sachs form factors. 
The point that concerns us here, the pion being a scalar 
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Fig. 2. Top panel: accepted [21] value of — 6F'(0), the would- 
be pion square charge radius, and lattice calculations [22] as 
a function of the pion mass. Bottom panel: the two terms 
of Eq. p4|) . that are way larger, and a cancelation must oc- 
cur to give the physical form factor derivative. It is not clear 
that one can extract the square charge radius, except for large 
quark masses where one could think of the lattice data in non- 
relativistic terms. 



target, is that the 1/m^ terms in Eq. for F'(0) are 
absent. 

The Breit (or brick-wall) frame is defined as that in 
which the target bounces off the virtual photon with op- 
posite incoming and outgoing momentum. As a conse- 
quence, EBroit = E^,.^it and the factor (p + p')72 (EE')^^^ 
in Eq. ([7]) and following becomes unity. 

Therefore in this frame. 



F(q^) 



(18) 



and therefore 



F(0) = (-f|f}ti°f^_ 



21 ' 



F'(0) = 2ik(-i 



F"(0) 



(-f|K4^^jOO_|-f) 



I 2 / ' 



(19) 



1 

6m^ 



F'(0). 
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In principle, the large term in F'(0) proportional to 
F(0) in Eq. ((T4| is absent, although the interpretation 
of the remaining in terms of a square charge radius is 
still not possible due to the interactions (Moeller opera- 
tors) and the necessary change of momentum (and hence 
frame) of the incoming and outgoing target. Moreover, the 
curvature in F"(0) is still affected by 1/mJ corrections. 
Therefore there is no clear advantage of interpretation of 
the lattice or experimental data in employing the Breit 
frame, at least in the case of the pion. 



(<P representing generically an integral over all fields). Since 
we are interested in exposing the explicit transverse glu- 
ons, the degrees of freedom of the intuitive Coulomb gauge, 
we divide the Faraday-Maxwell tensor in terms of chromo- 
electric and chromomagnetic fields 



E° = -9" A" - VA^'" + gf"'"=A*'A°'= 



V,A? - -gf 



j k 



(25) 



3 Poincare generators of QCD 



Coulomb gauge is fixed with the help of the Faddeev- 
Popov mechanism by first introducing a gauge-fixing term 
in the action 



Having motivated their phenomenological importance, in 
this section we derive the generators of the Poincare trans- 
formations in Coulomb gauge QCD. The method we follow 
is a modern path-integral approach in the first-order for- 
malism (with the path integration extending over phase 
space (q, p) [13], not only the field coordinates q), based 
on the work of Zwanziger [U and Watson and Reinhardt 
pS] . First we obtain the result for Yang-Mills theory in 
the functional approach. Then we derive the result for the 
boost generators appropriate for canonical quantization, 
in agreement with |26| . Finally we add the quarks to com- 
plete the Boost operators in canonically quantized QCD. 



2 = J V^exp {iSym + iSpp} , 
Sfp = / d^x [-C^ ■ A"" - c" V • D"''c''] . 



(26) 



3.1 The Hamiltonian and Momentum Generators of a 
pure Yang-Mills theory 

Yang-Mills theory is based on the renormalizable Lagrangian 
density of a spin-1, color-octet field A^(x), and written 
down in a gauge invariant manner in terms of the Maxwell 
tensor = 9^A^ - d^A^ + gP^'^A^A^ as 



Here c° are Grassmann (ghost) fields that keep track of 
the Faddeev-Popov determinant for now, involving the 
covariant-derivative operator D"** = 5"'' V + gf'''"^A'^ and 

is a Lagrange multiplier to assist with the Coulomb 
gauge constraint. 

As is known. Coulomb gauge fixing is not unique in 
non-Abelian gauge theories. Gribov [27| was the first to 
point out that there remain physically equivalent gauge 
configurations that are related by finite (as opposed to 
infinitesimal) gauge transformations. A further restriction 
must be imposed to the configuration space of gauge fields 
to a region where 



i7 = {A: V-A = 0| -V-D>0}. 



(27) 



5- 



YM 



d^x^YM with J^^ 



YM 



(20) 



The canonical stress-energy tensor conserved by transla- 
tional invariance of ^ym reads 



^YM 



(21) 



(22) 



This simple domain is still not totally free of Gribov copies. 
Instead one should consider the fundamental modular re- 
gion A [5^. It turns out though that functional integrals 
are dominated by configurations on the common bound- 
ary of fl and A [55] so that, in practice, it is enough to 
consider the domain defined by Eq. (P7)) . 

The canonical momentum associated to the transverse 
A" fields is introduced as an auxiliary field thanks to the 
following identity 



g'^'^being the Minkowski metric with signature H . 

Associated to space-time translational invariance are the 
conserved Hamiltonian and Momentum 



Ptt exp \ i \ d X 



exp <^ 
1 



d'^x-E'^ ■YI'S ^ 



TV -TV — TT • Hi 



(28) 



H 



YM 



d^x^YM and Vym 



YM- 



(23) 



In the path-integral framework the Green's functions 
are generated from the functional 



Z = y V$exp{iSYM} 



(24) 



Two more auxiliary variables f2, t will allow dividing the 
chromoelectric part of the action into the dynamical trans- 
verse potential of Coulomb gauge and a constrained lon- 
gitudinal part. One needs the further identity 

const = j VQ5 (V • TT + V^I?) 

= j V{n,T} exp|-i j d^xT'' (V • tt" + V2i7")|. 
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and a change of variables tt — > tt — Vi7 to arrive to the 
equivalent form of the Yang-Mills action 



->YM 



iB"-B°-r"V-7r" 
2 



(29) 



Wc now perform the integral over the ghost fields in 
Eq. (pS)) to recover the Jacobian of the change of vari- 
ables. Likewise we integrate over the Lagrange multipli- 
ers and r" to recover the delta-functions associated to 
the constraints, that enforce transversality of the physical 
fields (the Coulomb gauge condition V • A = V • tt = 0). 
Thereafter one can set to zero any terms involving these 
divergences in the action. 

Taking these details into account Z can be expressed 

as 

Z = j P^Det [- V • m^{yi - y)] 6 {V ■ A) (30) 
X(5 (V • tt) exp {iS} 



with 



S 



d^x 



-B" B" 



— TT • TT 

2 



2 



tt" • oi"A" + A°° (V • D"*"!?^ + gQl)] 



(31) 



and q'^ ~ ^adeji^d . color-chargc density carried by 

the gluons. 

The functional integral over AO can also be performed 
to yield the constraint equivalent to Poisson's equation in 
Quantum Electrodynamics in the form of a new delta- 
function on the path integral, 5 (-V • D'^^U^ - gp^). The 
Coulomb instantaneous potential of QED, l/|x — y|, in- 
verse of the Laplacian, is generalized in QCD to the inverse 
of the Faddeev-Popov operator 



be 



(32) 



and the constraint can be formally solved for il, now tak- 
ing the place of the scalar potential A^, 

= gM^^g^ (33) 

The associated delta-function factorizes in a useful way [H] 

5 (- V • D-^"/?^ - ge^) 

= Det [-V • D<54(x - y)] 6 {H'^ - gM'^'^gl) (34) 

to cancel the determinant in Eq. (j30p yielding 



Z 

where 



VAVtvS (V • A) 6 (V • tt) exp {iSo} 



So = I d*x 
1 



-Ib" -B" - ^tt" - tt" 
2 2 



(35) 

(36) 
(37) 



Clearly, the term tt" • A° is equivalent to pq in classical 
mechanics, and the remaining part can be identified as the 
classical Hamiltonian of pure Yang-Mills theory 



H 



YM 



yjrYM(x,t)d3 



J^YM = ^7r''(x,t) • 7rnx,t) + iB''(x,t) • B''(x,t) (38) 



+ -g2^,^M'«^(-V2)M-e^. 

We now take our expressions to the canonical func- 
tional quantization formalism. At this point tt ceases to be 
an integration variable and must be considered as the stan- 
dard classical canonical momentum tt" = 9^ym/9(A)°. 
After canonical quantization, the transverse field tt" has 
to be substituted by the momentum conjugate 11° to the 
transverse field A°. This is however not trivial, but the 
problem has already been solved in the past. 

Lideed, comparing this Hamiltonian with that of Christ 
and Lee [50] , or in recent papers [5T1I5^[551IM] , one notices 
the absence of the Faddeev-Popov determinant in our ex- 
pression. In order to introduce J^[A\ ~ Det(— V • D) we 
first note that the path integral representation in Eq. psp 
involves a Cartesian integration measure over gauge field. 
Let us consider, then, the corresponding quantized version 
of Eq. dig) 

J#YM = ^n'^^(x) • n''(x) + iB"(x) • b"(x) (39) 
+ ig2f'"''=A''(x) • n"t(x)M'"'(-V2)M''=|^ 

One is not entitled to use these operators (that we have 
covered with a tilde) as the canonical momentum since 

7^ — i(5/(5A°(x) (we repair this below in Eq. (^5]) ). In 
a Hilbert space where the scalar product of the complex- 
valued wavefunctionals !f'[A] is taken as "Cartesian": 



PAi^i [A]#2[A] . 



(40) 



(Note that in Eq. ^ we have fixed t = 0). 
We now, introduce 

tf^^A] = J-i/2[A].^,[A], 



(41) 



which is the wavefunctional one would find if, as Christ 
and Lee did, one had first quantized in A" = gauge 
and then transformed the gauge fields to Coulomb gauge, 
since in this variable change the scalar product picks up 
the determinant Eq. (|^IJ| 



{^M^ J VAJ[A]^*[A]^2[A] 



(42) 



that we would have instead absorbed in the wavefunc- 
tional. However, by doing the change of wave functional 
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Eq. (|41D one should at the same time preserve the expec- The quantized version in the flat gauge field configuration 
tation value of observables (!^i|(!)|#2) = (>f'i|0|!f2) with ^^^^^ 

(3 ^ ^-1/2^^1/2^ (43) ^YM = -J d3xn^-(x, t) VA^'^(x, t) (50) 

This has to be true, in particular for the Hamiltonian den- Finally, after applying the transformation 
sity. As a consequence 

Vym = J-'^^VymJ'^^ (51) 
and by considering Eq. (|45l) we obtain 



= ij-i/2ftat(x)n''(x)^i/2 + 1b"(x) • b"(x)(44) 
+ ig2j-^/^e^M'"'(-V2)M''=|^J^i/2. 



YM 



Y dW^(x,t)VA^^(x,t) (52) 



(in agreement with j26j and generalizing the momentum 
Replacing the operator n''(x) by the transformed one electrodynamics). 



n"(x) = J^^/2fl''(x)J-i/2 (45) 



3.2 Rotation and Boost Generators of pure Yang-Mills 



= n'^(x) + -^^^ theory 
2^A''(x) 

In this subsection we obtain the generators associated to 
we find that the quantum Hamiltonian density features the homogeneous Lorentz group. We again invoke Noether's 
the Faddeev-Popov determinant current for the field A'^(x), (with special transformation 

5A^ = l/2a."/^(^„^)>'^) 

^YM = -J-^il'^Wil^^) + -B (x) . B (x) (46) 

+ \iJ-'9\ji^'''{'-^'W^9%. ^^'^'^"^ ^ {5.A''^(^'^^)^x'^ (53) 



2' 



Now the dynamical operators in Eq. pS)) are A''(x) and 

Il'^(x) = -z5/(5A«(x), that do satisfy the equal-time com- ^^ere (^„^)f = b^g^i - b'-^g^i is the vectorial represen- 
mutation relation tation of the Lie algebra generator of S0(3, 1). Recalling 

Eq. we can write 

A^(x),77,^(y)l ^^b'^^{b,^-Y^s\^-y). (47) 



Y72 ]" ^ I- I Jym = £'^"'^ -I- S^"'' 



YM 



A prize to pay is that n°(x) is not a Hermitian opera- £^"5^ = x"5^]^ - x^^,^, (54) 
tor for the curvilinear scalar product. Indeed, n'^t(x) = 

^-in-(x) J ^ n'^(x) (see details in appendixEl) 6'^"^ = - ^Jf-A^-. 

We turn our attention to Vym defined in Eq. (^5)) . 

Taking into account Eq. ^ we can express Here is the orbital part whereas 6^'"" corresponds to 

the intrinsic spin. Note that Jym include both rotations 

/d'^xE^^fx t)VA^°fx t) (48) ^"^"^ boosts belonging to the Lorentz group. We will start 

y ^ ' by extracting the first. This is obtained from Eq. ([M]) by 
fixing /i = and taking the spatial part of the remaining 

The structure of the above expression does not difl^er from tensor. As a consequence the structure of the angular mo- 

those obtained for the case of a pure U(l)-gauge theory, mentum does not differ from an Abelian group except by 

New here is the sum running over the color indices. a sum running over the color indices 

We then split the chromoelectric field into transverse 

and longitudinal part E° = E^j. -|~ Vi?". From a classical Zym ~ J^ym + <5ym, 
point of view E^j. = — tt". Under these conditions 

PYM--/d3x{7r^'^(x,t)VA^'^(x,t) £YM = |d3x{xx [E^-(x,t)VA^^(x,t)]}, ^^^^ 

(49) 



+ a^i7'^(x,t)A'^°(x,t)} 



>YM 



J d^x{E^(x,t) X A''(x,t)} . 



Next, we integrate by part the second term of Eq. (|49)) and 
use the identity V • A° = 0. As a consequence, it vanishes where we have used = E"'. We again decompose 

identically and only the first term of Eq. (UHl) remains. E" — tt" + VJ7°. Terms arising from the longitudinal 
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part, VJ?° cancel each other. In the end, As for the longitudinal part, Eq. substitution of 

solutions of Ea. ((62)) and Eq. (|60| give 



2ym^- j d'xjx X [n^"(x,t)VA^"(x,t) 
- JT"(x,t) X A''(x,t)| . 



(56) i^ii = -g' j d3xV(f?^M'"')M-e^ . (65) 



After integration by parts = g^ / d^xgi^M'"' V(M''^£ip. 
Next we address the Boost generator in pure Yang- Therefore R\\ is a total derivative and vanishes identically. 
Mills theory, fixing ^ = a = in Eq. ([54]). This procedure Wrapping up, the combination of Eq. and Eq. 
leads to in Kym allows us to express it as 

JyM ~ X ^YM X ^YM +i> , {Oi ) ^ j ^3^^ J MtsQ T3a I 

and must be understood as the Hamiltonian density 



Kym = - / d-^x <j X ( -B° • B° + -tt" • tt" (66) 



in Eq. ([55)) . Quantizing again at t = 0, the corresponding _l_ l^g'']y[''a(_y2-j]y|ac^c \ j_ %.j^<^y[<^b 

charge is 2 ^ 2 " 



Kym — K-YM + -^ym, 



^ __ /"(^Sx/xfiB" B" + -tt" tt" To match the expression of Besting and Schiitte it suf- 

_/ " \ \2 2 ,p.„^ fices to add to Eq. a vanishing term proportional to 

+ .^''-(-V^lM-o^^ 1 ^"^ " = Ng^jd?^V{Q\M^-)M-^Ql, and choosing = 

2 8 / J -'^Z^' ^^^^ ^'^'^ corresponding boost operator for a fiat 

scalar product can be written as 

i^YM=/d3xE''(x)Ag(x). ^ [ { fl 1" ^ \ 

Returning to Eq. the classical Euler-Lagrange equa- ^™ ^ j i ' ^ ^ 2^"^ ' J " ''^''^ 

tion of motion associated to Ag, Gauss's law, is a con- 2 

straint imposed on the quantization (above in Eq. ([M]) we - —Q'^^A^'^{^|^:i^^^)W■''Ql + -n°'''M°^gg -|- 

used it in the form of the classical Poisson equation) ^ ^ ^ 

^ . u . + |ea'^M''''n^]. . (68) 

V • - gf'^''"A'' • E'^ = D'''' • E'' = (59) 2 ^ J 

Substitution of E" = - A - D'^^'Aj," in the first term of Eq. Making the change Kym = >7"^/^Kym J^^^^ and mind- 

(|59p returns ing the correct ordering of the various, non-commuting 

— V • D'^'^Aq ~ ggg (60) operators we obtain 

where we have used the trasversality condition (V • A"^ = , /"^fli- - 1-.- 

0) and identified = f'''"=A'' • tt". Once more we de- Kym = - j d'xj-J" H^xH^ -f- -B^xB^ 
compose the electric field into transverse and longitudinal 

P""'*' E'^ = -7r'^ + Vr2'^ (61) -\g^J-^glJM^'^{du^d,)M<^^gl (69) 

in the middle term of Eq. ([SS]) to derive the equation _^ -gJ'^^JI^ JM^^ g° + -gJ~^g° J'M^''ri'^| 
D"'' • Vl?^ = g(?° (62) 



Eq. (j69p coincides with the expression derived in |26j . To 
(at this point Aq and Q are interchangeable, but we keep derive the above expression we have replaced the operator 
the distinction for another instant). Taking into account fi'^iy^) by the transformed one (see Eq. (03). 



Eq. ((6T|) the spin piece of the boost, Eq. ([58|. becomes 

i^YM = il± + i^ii , 3 3 Inclusion of the Quark sector 

.^j^ = — g / d'^x7r"M"^£i^, To complete the discussion on Quantum Chromodynamics 

^ ' we include the matter sector interacting with the gauge 
field. This is achieved by the matter Lagragian 

=2'Matter = q'(«7^9M-m)q^ + g^V^A^q^ (70) 
We express ^± in a more symmetric form ^ 

with q = q^'7° and A"/2 representing the SU(3) Gell-Mann 



^W^eJ d3x(Vr2'')A°. 



= -f / d^xTT'^M''^,^, - I J d^g^W^n^ (64) ^,'^trices. 
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Its contribution to the stress-energy tensor is 



"^Matter 



(71) 



i^'-1^d^i( - g^'^Wi-i^dx - m)q^ - (72) 
Y 



The conserved charges associated to space-time transla- 
tional invariance are then 



HMattcr = / d^X^iOOtt,, = / d^X {^^Dirac + ^nt } , 



(73) 



PMatter = / d^x^^!,,,,, = / d3xq^t(x) (-^V) q^(x) 
where J^imc is the free Dirac Hamiltonian density: 

^Dirac(x) = q^t(x)(-ia • V + /?m)q^(x), (74) 



whereas 



^„t(x) = gq^^x) ya • A'^(x)q^(x) (75) 



A" 



arises due to the interaction between the gauge field and 
Quarks. The last term of the above equation introduces 
the quark color charge density g'^ = q^^'^q^ that must be 
summed to the right hand side (source) of Poisson's equa- 
tion in Eq. (p3|) and following. As a consequence ^• 
^ M'"''^'' with = + ^q- The second term of 
Eq. ((75|) adds up to the second term of Eq. (|46| and both 
together have the same functional form as the latter but 
with replaced g°' = g"^ + g^- 

The complete QCD Hamiltonian and momentum den- 
sities are the sum of both 



HqCD = HyM + HMatter, 
"PqCD = T-'yM + "^Matter- 



(76) 
(77) 



where Hym is the corresponding operational version of 
Eq. dSHD with ,^YM given by Eq. (gS]) In turn, Eq. ((771) 
combines Eq. ([52]) and Eq. (|73)) . 

We now turn to the angular momentum and boost op- 
erators. The conserved current due to quarks is obtained 
from the tensor 



"-'Matter ~ ^ "^Matter 



^Matter + 4 q 



* -i aB f^-SfMattor 
-q cr . 



(78) 



Here 5ifMattcr/c'(9^qO = ^^7^ whereas 9^Mattor/9(5^q^ 
—ij^q therefore 



•^MattorV-^i ~ "^Matter "^Matter' 



■^Matter 

■ 



6 

^ ^Matter ^ ^AA.t 



Matter' 



(79) 



S 



i^^{7^a"/5}q^. 



By fixing fi ~ and taking the remaining indices as 
spatial, we obtain the quark angular momentum, ^Matter = 
l/2e*^'= / d3x(x.'-^j0^^,,„ - iq^t^.fcq^), ^here Tgf,,,,, is the 
momentum density in Eq. (|73p . Therefore 



al Matter — ^ 



Matter 



'Matter: 



'Matter 



y'd3x{q^t(x)[xx(-»V)]q^(x)}, (80) 



•Matter 



/dW^(x)(ii;) q^(x), 



with Ui ~ \f-ijk(^jk- 

Following a procedure similar to that used in deriv- 
ing Eq. ((57)) one obtains the Boost density easily. The 
generator is obtained integrating over the (t = 0) surface 



Matter ~ / '^^^ { "Matter 2 

.f^Matter = ^Matter + -^Matteri 

/Cwatter = - j d^xx {-ict ■ V + /?m) q^(x) 

+ |q^t(x)A''a- A'^(x)q^(x)}, 



iq£tcrOiq<'} . Therefore 



(81) 



Matt cr 



/dWt(x)(f«) q^(x) 



with ai = icFoi- For complete QCD, one needs to substitute 
again g^ by g"' = gg + g^ in Eq. ((5S|) , the boost generator 
for pure Yang-Mills theory. Finally one finds 



5qCD — 5yM + -^Matter: 
KqcD ~ Kym + ^^Matter- 



(82) 
(83) 



(keeping in mind that gg~j'g°' = gg + gq)- 

We also record the (possibly Bogoliubov rotated) nor- 
mal mode expansions of the dynamical fields 

A''(x) = / [a''(k) +a'^t(_k)l e'k-x^ 

J v2wk 



n (x 



dky^[a''(k)-a'^t(_k)] 
q^(x)=^ /dke"^- (^<AUkA4 

q^t(x) ^Y.[ ^ke"'" '' (<K + i?i,,Vl 



kA ^ -kA 



(84) 



kA 



with dk = d^k/(27r)'^ and Wk — |k|. In momentum space 
Eq. (|37]) reduces to 



a^(k),a^^^(k') ^{27TrS'''[S., 



^^)5^(k-k') (85) 



Transversality translates into k • a" (k) = k • a"t (k) . The 
operators a" (k) and a"''' (k) include the gluon polarization 
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a" = J2r=i 2 Ci-(k)a°(k) and the usual creation a°t(k) and 
annihilation a° (k) operators (the normalization is as usual 
eiel = 6rs, el.e>. = S'^ - k^k^k^). 

Similarly, quark creation and annihilation operators 
satisfy the anticommutation relations 

(86) 

(others zero). 

The renormalization of the boost generators also de- 
serves a comment. It is clear that radiative corrections to 
the classical, conformal theory are infinite and the quan- 
tum theory is only well defined after a scale is chosen. 
Regularization and renormalization in fixed-gauge Hamil- 
tonian dynamics is arduous and is being pursued indepen- 
dently with attention to the Slavnov- Taylor identities |35j . 
One should like to find a Lorentz-invariant regulator that 
would leave the Poincarc algebra intact (which could be 
established by looking into Schwinger's condition). What- 
ever strategy one adopts to regularize the Hamiltonian, for 
example a lattice regularization, one can similarly apply 
it to the boost operator. A large part of the same, com- 
ing from x'H, is automatically regulated and requires no 
additional counterterms other than those in H, as the ad- 
ditional power of x' lowers the mass-dimension. This leaves 
the spin terms like .ft, and in principle, one should expect 
mass counterterms of equal or lower mass dimension than 
those spin terms to appear. The issue is postponed to fu- 
ture work. 



4 (/)-radiative decays 

4.1 The boosted wavefunctions of decay products 

We consider again the radiative decay -> 777 as an exam- 
ple, but it is obvious that the discussion is general. The 
boosted meson with velocity v is given (in terms of the 
rapidity ( defined in Eq. ^) by 



(87) 



(K • 



We want to establish that a gluonium term, not taken into 
account in past analysis, arises just because of the change 
of reference frame, even if \r]o) contained only purely qq 
configurations. Therefore we need to see what pieces of 
the QCD boost operator connect the qq and gg Fock sub- 
spaces. 

At linear order in v, the two spaces arc indeed discon- 
nected, since the interacting part of K that could effect 
the change is proportional to the color charge densities in 

quarks and gluons respectively, QqQg oc l/2q'''A''qf'''"^A • 

^ c 

II , and upon closing the quark line, color arithmetic sets 
the contribution to zero since Tr(A) ~ 0. One needs then 
to resource to the second order term, K'^, where the square 



of the instantaneous Coulomb potential contributes, as 
does the square of the gauge coupling to transverse gluons. 
We focuse on the latter 



d3xzq^1'(x)A'^a- A (x)q^(x) 



to show the phenomenon, not because the contribution 
of the a ■ A Hamiltonian is larger than the gg piece, but 
just because of simplicity: this contribution is the one with 
the smallest number of loops that turns out not to van- 
ish, which we here demonstrate. An added bonus is that, 
although formally in Coulomb gauge, the piece we cal- 
culate will be present in other gauges since it involves 
only the transverse gluons. The calculational methods are 
analogous to those of the estimate of the glueball width 
in Coulomb-gauge models of QCD }36j . 

With V ~ 0.55 for our particular case, ^ ~ 0.15. For 
the rest-frame state is sufficient to take the color singlet, 
fiavor singlet (through mixing) component of the ry meson 



1^0; 



= sin(0p)y"d3k Y°(f2k)^„„(|k|) (89) 
X E(^/^i^-A^2| 00^ (-1)1/2+'^^ <,^lk,jO) 



with color wavefunction S^^^^ /{NcY^^ , and flavor wave- 
function understood to be the singlet projection. 

Substituting the normal-mode expansion Eq. (|M|) for 

^ a 

the fields q and A in Eq. (|88|) . one obtains, in terms of 
BCS spinors U, V, a gluonium component of the 77 moving 
with velocity v, 



\Vv) = g^y smi^ipp 



dq y dp y d^k J d\' 

dp^ (3) (k + p _ q)) (^^(3) (k' + q - p) 



(90) 



1 1 

—jll — 

An 2^WkWk' \ 2 2 



M2 



00^ (-l)l/2+M2 



The color algebra is easy since Tr(A°A'') = 2(5"''. The spin 
work however is a little more tedious. Let us adopt a def- 
inition of the spinors in terms of a BCS angle, defined as 
Sk EE sin^(|k|) = with E(|k|) = [k2 + m(|k|)2]i/2 

and Ck = cos(/3(|k|) = p^TTkTT- spinors read then 



E(|k|)- 



UkA 



V_kA = 



21/2 



21/2 



(1 + Sk)l/2XA 



(l-5k 



\l/2 



or ■ nkXA 



(1 - Sk)^/^cr • nkXA 



(91) 
(92) 



in terms of (bidimensional) Pauli matrices and Pauli spinors 
satisfying the closure relation J^x XxXx ~ l2x2- Here nk = 
k/|k| is a unit vector. 



Gomez Rocha, Llanes-Estrada, Schiitte and Villalba: Coulomb-QCD boosts and applications 



11 



The spin combination needed yields 



tions being integrated. To proeeed we write the overlap as 



A 

A 

+ (1 - Sq)cpcr • Uqcr^cr ■ HpX^iXj^^cr^ cr • Hq- 



(1 + Sq)CpCr-'XMlXL'^ ■ "P'^'' - 
(1 + Sp)CciCr^ Xt^iXU'^^ cr • nqj XA 



(93) 



Then we take into account that the Clebsch-Gordan 
coefficient in Eq. (|M|) is 



(^Mi ^-M2| 00^ (-1)1/2+^^ = -i=5,,,. 



(94) 



which forces a trace over the Pauli matrices. Employing 
cyclicity of the trace and Ti{a^a^ a'') = 2ie^^^, and re- 
turning to the running mass m(|p|) and energy E(|p|), we 
have 



\Vv) = y(47ras)sin(0p)ydq ^dp 

1 1 q'm(|p|) -p*m(|q|) 



(95) 



4^ V2 E(|q|)E(|p| 



X da 



^ ,tai \ (_p, \ I 1 Jaj 



<4 (-^qJ 



q-p 



<4 lo)- 



p-q 



where as = g^/(47r) is the strong fine structure constant. 



4.2 The gluonium content 



In calculating matrix elements of the boost operator one 
encounters unusual features. The presence of turns, in 
the usual momentum representation, in a derivative re- 
spect to fc', that affects in the end the (5-functions from 
the field (anti)commutators. While more than one such 
momentum-conservation functions appear, the following 
construction can be found often: / dkF{k)5{k)d\5{k) that 
we at present ignore how to handle. If we ignore the term 
with second derivatives acting on the Dirac delta func- 
tions, we can complete the rest of the calculation by the 
usual per-parts integration method that systematically 
transfers derivatives from the (5-functions to the other func- 



{r^l\r^l) = {Airasf ( ^ ) sin^<pp; 



X (m(lpl)q^ - m(|q|)p')(m(|p'|)q ^' - m(|q'|)p ^ 
d^kid-^k^d^k^d^k^, ^ (3) 

— ^—TJ^9^1 S^'^ (ki + q - P) 

X ak|<5(3)(k2 + p - q)d^>^j^^\K + q' - p') 
X \','S^'H^2 + p' - q')(27r)« {<5(^)(ki - k'l) 

X t^'^'(ki)<5(3)(k2 - k^)t"'(k2) + <^^''(ki - k^) 
X t^''(ki)^(3)(k,_k'i)t'J"'(k2)}. 

There are several things to remark concerning this ex- 
pression. First, note that t'-'(k) = {S^^ — k'k^/k^) is the 
transverse projector in momentum space which ensures 
transvcrsality. Since the suite of Kronecker deltas will force 
ki = — k2, and due to the spin contraction with anti- 
symmetric e*^' symbols, the terms quartic in fc„ will not 
contribute. It is also convenient to define a shorthand no- 
tation regrouping the meson wavefunction and the quark 
propagator pieces, 

G(p,q).^.„(|q|)^^^»-^^ (97) 

(where the second wavefunction depending on |p| arises 
due to the second contributing spin combination). 
Finally, as explained above, 

5(3) (ki + q - p)ak3 5 (ki + q - p) X [. . . ] 
^ ,5(3) (ki + q - p)5(3) + q _ p)dp. [. . . ] 

ignoring the effect of the derivatives over other delta func- 
tions, that would normally lead to the SW'^S term. Due to 
invariance under rotations, one can also substitute deriva- 
tives by the Laplacian d'^/d^ 

After isolating the (5's, one encounters two factors of 
the type 



1 



1 



9 "1 [2wki2a;k; 
that are handled by means of 

1 



-t^' (ki)t^' '(k2) 



t'^(k) 




-I 


\ 4cjk 2 


W|k|2 \ 





Since there are two insertions of the boost operator, this 
is a three-loop calculation. The integration variables can 
be chosen as q, p and q'. The fourth, p' becomes a linear 
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combination. There are two terms with a different value 
of p' in each, that we will denote with a subindex 

p'l = p q + q' 
Pa = p + q + q' . 

We can then write down our partial calculation for the 
overlap as 



(0|0) 



dq / dp 



sm 







1 1 



(99) 



with spin work that can be obtained by (.^q_p 

(q-p)/|q-p|) 



^ijl^i'j'l' 



{S"' -3Z\M^')(5"' -3Z\'Z\'')G"'G''' ^ (100) 



-2G-G' + 3(A-G)(A-G') 

and similar relations, to yield (uj = Wq_p) 

2 



LxP- \ 4a; 



|q-p| 



[G(p,q).Ziq_^^01) 




:,q).G(p;,q') + 2[G(p,q).Aq_p] 



|q-p| 



(-G(p, 



[G(pl,q')-Aq_p]) + (-2G(p,q> 
G(p;,q') + 3(G(p,q)-Aq_p 



-pr 

_p)(G(p;,q')-Aq-p)} 



and with obtainable replacing the argument of G(pi , q') 
by the second solution to the momentum 5's G(p2,q'). 
With all these factors, one can proceed to a numerical 
evaluation of Eq. (|99p . Given the now explicit invariance 
under rotations of the integrand, the nine-dimensional in- 
tegral can be reduced by choosing the z axis along p, and 
the X axis so that p, q are in the xz plane. The remain- 
ing integration variables are then the three moduli |p|, 
|q|, |q'|, and the three integration angles 0q, l?q'. This 
six-dimensional integral can be carried out with the help 
of VEGAS, the standard adaptive Montecarlo algorithm 
|37| . The running quark masses m(|q|) and gluon (disper- 
sive) energies w(|k|) appearing satisfy their respective gap 
equations. We take them to be solutions of the truncated 
versions in [35] and [3^, as later work, especially in the 
gluon sector, has improved the analysis but not changed 
the numerical results qualitatively. 

Exploration of the six-dimensional numerical integral 
9, defined by 



(r;8k«) = (47ra.)^(^ 



sm 



(102) 



and taking into account conceivable uncertainties in the 
quarkonium 77 wavefunction, quark mass gap function, and 
gluon dispersive function, allows us to make an order of 
magnitude estimate (note 3, after extraction of the vol- 
ume factor (27r)^5'^'^^(0) is a pure, dimensionless number) 
log 3 = 2.6^2 0- ^^'^ prefactor shows the explicit depen- 
dence in the pseudoscalar mixing angle (through which 
the singlet component of the t] arises), which we take at 
39 degrees, the rapidity C — 0.62 of the boosted t] meson, 
and as at some low scale. This last value is not known 
with certainty, but since it is estimated to be 0.4 in t de- 
cays, we obtain a lower bound to the prefactor of 3 which 
is 0.4. 

Putting all together, we obtain 



3+0.5 
'-2.0 



(103) 



which means that (?7^|f?5) is larger than 1 by somewhat 
more than one sigma. This is of course not possible, as 
the sum of all coefHcients of the wavefunction in Eq.Q 
5^- c? = 1, but one has to take into account that this is 
only a partial calculation including only one of the relevant 
terms of the boost operator, and that model uncertainties 
creep in the choice of wave and gap functions. 

What one should conclude is that there is no reason 
why {T]%\r]^) should be orders of magnitude smaller than 
one. By the size of the controled factors one could guess 
it to be of order 0.1. That is, we believe we can argue 
with certain confidence that a gluonium component of the 
T] meson is induced by boosting it to a moving reference 
frame, as in (\) radiative decays, even if this component 
was not present in the rest-frame wavefunction. 



5 Summary 

In this work we have derived the boost generators of Quan- 
tum Chromodynamics quantized in Coulomb gauge in the 
path integral formalism. These confirm the result of Best- 
ing and Scliiitte for pure Yang-Mills theory, to which we 
add the fermion fields to complete the boost operator of 
canonically quantized QCD in Coulomb gauge 

Kqcd = - y" d^x I i J-in» JxH'' + ie'^xB'' (104) 

+ ig^-in^^M"''^'' + ji^'^ii'' 

- J d^xx {q^t(x) (^ia . V + /3niq) q^(x) 
+ gq^t(x)^«.A'^(x)q^(x) 



dWt(,) ' 



We have wanted to call the attention of the commu- 
nity to the dynamical nature of these operators. We have 
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focused on two aspects of current interest for the interpre- 
tation of hadron data in terms of their quark and gluon 
constituents. The first is the lack of a charge-distribution 
interpretation of the pion form factor. The operator ma- 
trix elements that one obtains from such experiments are a 
much more complicated generalization. The second is the 
fact that Fock-space expansions of hadron wavefunctions 
are tied to the rest frame of the hadron, and they change 
with its velocity. This is of current interest for the theo- 
retical understanding of pseudoscalar mesons produced in 
(f) radiative decays at Frascatti. 

Many more applications come to mind where theoreti- 
cal interpretation is complicated by the dynamical nature 
of the boost operators. However there is no known prac- 
tical way of handling them to obtain even model results. 
We look forward to any progress in this direction by lattice 
methods or other means, n 



The introduction of a unit factor 1 = J' J' ^ in front of 
the bracket of the above expression leads to 



VAJ 



(107) 



It is to notice that 

(lfl|fl''t(x)|*'2) = (n-(x)lfi|tf'2> 



(108) 



= / -DAJ 



and we compare Eq. (|107p with Eq. pOSp to find that with 
the curvilinear scalar product, the canonical momentum 
is not Hermitian, n'''l'(x) = J~^Il''(x.)J . 
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7. 



A The Fadeev-Popov Determinant and the 
Hermiticity of n"(x). 

Let us first recheek for the reader not familiar with func- 
tional quantization how Hermiticity of the canonical mo- 
mentum comes about, in total analogy with elementary 
quantum mechanics. 



(tf'i|ri'^(x)|iZ'2) = / VAJ>F*[A] 



1 



(5A«(x) 



IP's [A] . (105) 



Again integrating by parts 

{^i\n°'{x)\W2) ^ - [vA 



1 s 

i (5A°(x) 



MM 



VA 



n"(x):7<fi[A] MM ■ (106) 



9. 
10. 
11. 
12. 
13. 
14. 

15. 
16. 
17. 



^ Of course, Light Front quantization or explicitly covariant 18. 
approaches such as the Dyson-Schwinger equations are free of 
this problem by construction but they have other issues of their 19. 
own that need to be overcome. Coulomb gauge QCD remains 20. 
an attractive avenue of thought for non-perturbative problems 
since the physical Hilbert space makes transparent the appli- 21. 
cation of the variational principle. 
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